I Claim: 

1 LA method for determining the state of stability of an electrical grid having n nodes, 

2 comprising the steps of: 

3 a. embedding load flow equations (L) representing the electrical grid in a parametric 

4 homotopy (L(s)) that goes continuously from a 0-case (L(0)), in which all voltages are equal to 

5 the nominal voltage and there is no energy flow in links of the grid, to an objective case (L(l)) 

6 representative of the grid in the condition for which stability is to be determined; 

7 b. developing in power series values of the load flow equations' unknowns in the 

8 parameters of the parametric homotopy (L(s)) in a neighborhood of the 0-case value of each 

9 parameter; 

10 c. computing a continued fraction approximation to the power series coefficients 

11 produced in step b; 

12 d. evaluating the n-order approximant of the continued fraction approximation produced 

13 in step c for the power series coefficients produced in step b to provide a solution to the load 

14 flow equations (L); and 

15 e. displaying the solution to the load flow equations as a measure of the state of stability 

16 of the electrical grid. 

1 2. The method of claim 1, further comprising the steps of: 

2 prior to said embedding step, receiving data from a supervisory and data acquisition 

3 system representative of conditions of the electrical grid, and forming said load flow equations 

4 (L) from said data. 
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1 3. The method of claim 2, further comprising the steps of repeating said receiving step 

2 and steps a through e continuously to provide a continuous, real time estimation of the stability 

3 of the electrical grid. 

1 4. The method of claim 3, further comprising the steps of confirming that a set of 

2 voltages and flows contained in said solution to said load flow equations (L) are representative of 

3 a physical electrical state. 

1 5. A method of measuring load flow in a power generating system having an electrical 

2 grid comprised of n nodes, comprising the steps of: 

3 a. embedding load flow equations (L) representing the electrical grid in a parametric 



4 homotopy (L(s)) that goes continuously from a 0-case (L(0)), in which all voltages are equal to 

5 the nominal voltage and there is no energy flow in links of the grid, to an objective case (L(l)) 

6 representative of the grid in the condition for which stability is to be determined; 



7 b. developing in power series values of the load flow equations' unknowns in the 

8 parameters of the parametric homotopy (L(s)) in a neighborhood of the 0-case value of each 

9 parameter; 

10 c. computing a continued fraction approximation to the power series coefficients 

11 produced in step b; 

12 d. evaluating the n-order approximant of the continued fraction approximation produced 

13 in step c for the power series coefficients produced in step b to provide a solution to the load 

14 flow equations (L); and 
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e. displaying the solution to the load flow equations as a measure of the load flow in the 
power generating system. 



1 6. The method of claim 5, further comprising the steps of: 

2 prior to said embedding step, receiving data from a supervisory and data acquisition 

3 system representative of conditions of the electrical grid, and forming said load flow equations 

4 (L) from said data. 



1 7. The method of claim 6, further comprising the steps of repeating said receiving step 

2 and steps a through e continuously to provide a continuous, real time measure of the load flow in 

3 the power generating system. 



1 8. A method of measuring load flow in a power generating system having an electrical 

2 grid, comprising the steps of: 

3 a. generating a mathematical model of a known, physical solution to the load flow 

4 equations (L) in which all voltages are equal to the nominal voltage and there is no energy flow 

5 in links of the grid; 

6 b. using analytical continuation to develop a mathematical model of the current, physical 

7 solution to the load flow equations representing the current load flow in the power generating 

8 system; and 

9 c. displaying the physical solution to the load flow equations as a measure of the load 
10 flow in the power generating system. 
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9. The method of claim 8, said generating step further comprising developing a power 
series expansion of all quantities in a parametric homotopy (L(s)) formed from said load flow 
equations (L) in a neighborhood of the 0-case value of each quantity. 



1 10. The method of claim 9, further comprising using algebraic approximants to 

2 determine the sum of all coefficients of said power series for the load flow equations 

3 representative of the physical current load flow that is to be determined. 

1 1 1 . A system for measuring load flow in a power generating system having an electrical 

2 grid, said system comprising: 

3 a supervisory control and data acquisition system adapted to collect data from said 



4 electrical grid indicative of electrical conditions in said electrical grid, said supervisory control 

5 and data acquisition system being in communication with a microprocessor-controlled energy 

6 management system, said energy management system further comprising executable computer 

7 instructions to: 



8 a. process said data received from said supervisory control and data acquisition 

9 system into load flow equations (L) representing the electrical grid; 

10 b. embed said load flow equations (L) in a parametric homotopy (L(s)) that goes 

11 continuously from a 0-case (L(0)), in which all voltage are equal to the nominal voltage and 

12 there is no energy flow in links of the grid, to an objective case (L(l)) representative of the grid 

13 in the condition for which stability is to be determined; 



14 c. develop in power series values of the load flow equations' unknowns in the 

15 parameters of the parametric homotopy (L(s)) in a neighborhood of the 0-case value of each 

16 parameter; 

17 d. compute a continued fraction approximation to the power series coefficients 

18 produced in step c; 

19 e. evaluate the n-order approximant of the continued fraction approximation 

20 produced in step d for the power series coefficients produced in step c to provide a solution to the 

21 load flow equations (L); and 

22 f. display the solution to the load flow equations as a measure of the state of 

23 stability of the electrical grid. 



